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Planar quadratic differential systems occur in many areas of applied mathematics. Although 
more than one thousand papers have been written on these systems, a complete understanding 
of this family is still missing. Classical problems, and in particular, Hilbert's 16th problem 
[Hilbert, 1900[ [Hilbert, 1902 , are still open for this family. In this article we make a global 



study of the family QTN of all real quadratic polynomial differential systems which have a semi- 
elemental triple node (triple node with exactly one zero eigenvalue). This family modulo the 
action of the affine group and time homotheties is three-dimensional and we give its bifurcation 
diagram with respect to a normal form, in the three-dimensional real space of the parameters 
of this form. This bifurcation diagram yields 28 phase portraits for systems in QTN counting 
phase portraits with and without limit cycles. Algebraic invariants are used to construct the 
bifurcation set. The phase portraits are represented on the Poincare disk. The bifurcation set is 
not only algebraic due to the presence of a surface found numerically. All points in this surface 
correspond to connections of separatrices. 
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1. Introduction, brief review of the litera- 
ture and statement of results 

In this paper we call quadratic differential systems 
or simply quadratic systems, differential systems of 
the form 

X = p{x,y), 

y = q{x,y), 



(1) 



where p and q are polynomials over M in x and y 
such that the max(deg(|?),deg((?)) = 2. To such 
a system one can always associate the quadratic 
vector field 

d d 

X=p— + q—, (2) 
ox oy 

as well as the differential equation 



ix — pdy = 0. 



(3) 



The class of all quadratic differential systems (or 
quadratic vector fields) will be denoted by QS. 
We can also write system ([T|) as 



V 



Po+Pi{x,y) +P2{x,y) =p{x,y), 
-- qo + qi(.x,y) + q2{x,y) = q{x,y), 



(4) 



where pi and qi are homogeneous polynomials of 
degree i in (x, y) with real coefficients with p^+^i 
0. 

The complete characterization of the phase 
portraits for real planar quadratic vector fields is 
not known and attempting to topologically classify 
these systems, which occur rather often in applica- 
tions, is quite a complex task. This family of sys- 
tems depends on twelve parameters, but due to the 
action of the group G of real affine transformations 
and time homotheties, the class ultimately depends 
on five parameters, but this is still a large number. 

The goal of this article is to study the class 
QTN of all quadratic systems possessing a semi- 
elemental triple node. By a semi-elemental point 
we understand a singular point with zero deter- 
minant of its Jacobian, but only one eigenvalue 
zero. These points are known in classical litera- 
ture as semi-elementary, but we use the term semi- 
elemental introduced in [Artes et al., 2012] as part 
of more consistent definitions associated to singular 
points, their multiplicities and Jacobians. 

The condition of having a semi-elemental triple 
node of all the systems in QTN implies that these 
systems may have another finite point or not. 



For a general framework of study of the 
class of all quadratic differential systems we re- 
fer to the article of Roussarie and Schlomiuk 
[Roussarie fc: Schlomiuk, 2002| . 

In this study we follow the pattern set out in 
[Artes et al., 2006] . As much as possible we shall 
try to avoid repeating technical sections which are 
the same for both papers, referring to the paper 
mentioned just above, for more complete informa- 
tion. 

In this article we give a partition of the class 
QTN into 63 parts: 17 three-dimensional ones, 29 
two-dimensional ones, 15 one-dimensional ones and 
2 points. This partition is obtained by considering 
all the bifurcation surfaces of singularities and one 
related to connections of separatrices, modulo "is- 
lands" (see Sec. 16. 3p . 



A graphic as defined in Dumortier et al., 1994 



is formed by a finite sequence of points ri, r2, . . . , r„ 
(with possible repetitions) and non-trivial connect- 
ing orbits 7i for i = 1, . . . ,n such that 7^ has r j 
as a-limit set and rj+i as a;-limit set for i < n 
and 7„ has r„ as a-limit set and ri as w-limit set. 
Also normal orientations nj of the non-trivial or- 
bits must be coherent in the sense that if Jj-i has 
left-hand orientation then so does jj. A polycycle 
is a graphic which has a Poincare return map. For 
more details, see [Dumortier et al, 1994[ . 

Theorem 1.1. There exist 28 distinct phase por- 
traits for the quadratic vector fields having a semi- 
elemental triple node. All these phase portraits are 
shown in Fig. \^ Moreover, the following state- 
ments hold: 

(a) There exist three phase portraits with limit cy- 
cles, and they are in the regions Vq, Vi^ and 
555; 

(b ) There exist four phase portraits with graphics, 
and they are in the regions 5S4, 7 Si, I.3L2 and 
5.7Li. 
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Fig. 1. Phase portraits for quadratic vector fields with a semi-elemental triple node. 
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From the 28 different phase portraits, 9 occur 
in 3-dimensional parts, 13 in 2-dimensional parts, 
5 in 1-dimensional parts and 1 occur in a single 
0-dimensional part. 

In Fig. [1] we have denoted with a little disk the 
elemental singular points and with a little trian- 
gle the semi-elemental triple node. We have plot- 
ted with wide curves the separatrices and we have 
added some thinner orbits to avoid confusion in 
some required cases. 

Remark 1.2. The phase portraits are labeled ac- 
cording to the parts of the bifurcation diagram 
where they occur. These labels could be dif- 
ferent for two topologically equivalent phase por- 
traits occurring in distinct parts. Some of the 
phase portraits in 3-dimensional parts also occur 
in some 2-dimensional parts bordering these 3- 
dimensional parts. An example occurs when a 
node turns into a focus. An analogous situation 
happens for phase portraits in 2-dimensional (re- 
spectively, 1-dimensional) parts, coinciding with a 
phase portrait on 1-dimensional (respectively, 0- 
dimensional) part situated on the border of it. 

The work is organized as follows. In Sec. [5] we 
describe the normal form for the family of systems 
having a semi-elemental triple node. 

For the study of real planar polynomial vector 
fields two compactifications are used. In Sec. [3] we 
describe very briefly the Poincare compactification 
on the 2-dimensional sphere. 

In Sec. [Hwe list some very basic properties of 
general quadratic systems needed in this study. 

In Sec. we mention some algebraic 
and geometric concepts that were introduced in 
ISchlomiuk et al, 2001^ |Lhbre et al, 2004| involv- 
ing intersection numbers, zero-cycles, divisors, and 
T-comitants and invariants for quadratic systems 
as used by the Sibirskii school. We refer the reader 
directly to [Artes et al., 2006] where these concepts 
are widely explained. 

In Sec. El using algebraic invariants and T- 
comitants, we construct the bifurcation surfaces for 
the class QTN. 

In Sec. [7] we introduce a global invariant de- 
noted by X, which classifies completely, up to topo- 
logical equivalence, the phase portraits we have ob- 
tained for the systems in the class QTN. Theorem 



17.61 shows clearly that they are uniquely determined 
(up to topological equivalence) by the values of the 
invariant I. 

2. Quadratic vector fields with a semi 
elemental triple node 

A singular point r of a planar vector field X in 
is semi-elemental if the determinant of the matrix 
of its linear part, DX{r), is zero, but its trace is 
different from zero. 

The following result characterizes the local 
phase portrait at a semi-elemental singular point. 



Proposition 2.1. lAndronov et al., 1973\ 
Dumortier et al., 2000^ Let r = (0,0) be an 
isolated singular point of the vector field X given 
by 

X = A{x,y), 
y = Xy + B{x,y), 



(5) 



where A and B are analytic in a neighborhood of 
the origin starting with at least degree 2 in the 
variables x and y. Let y = f{x) be the solu- 
tion of the equation Xy + B{x,y) = in a neigh- 
borhood of the point r = (0,0), and suppose that 
the function g{x) = A{x,f{x)) has the expression 
g{x) = ax" + 0(3;"), where a > 2 and a ^ 0. So, 
when a is odd, then r = (0, 0) is either an unsta- 
ble multiple node, or a multiple saddle, depending 
if a > 0, or a < 0, respectively. In the case of 
the multiple saddle, the separatrices are tangent to 
the X-axis. If a is even, the r = (0,0) is a multiple 
saddle-node, i.e. the singular point is formed by the 
union of two hyperbolic sectors with one parabolic 
sector. The stable separatrix is tangent to the posi- 
tive (respectively, negative) x-axis at r = (0,0) ac- 
cording to a < (respectively, a > Oj. The two 
unstable separatrices are tangent to the y-axis at 
r = (0,0). 

In the particular case where A and B are real 
quadratic polynomials in the variables x and y, a 
quadratic system with a semi-elemental singular 
point at the origin can always be written into the 
form 

X = gx"^ + 2hxy + ky"^, 
y = y + + 2mxy + ny^. 



(6) 



By Proposition 12.11 if y 7^ 0, then we have a 
double saddle-node sn(2) , using the notation intro- 
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duced in [Artes et al, 2012| . Otherwise, if g = 
and £ ^ 0, then, if i < 0, we have a triple node n(3) 
and, if £ > 0, we have a triple saddle S(3). The mul- 
tiplicity of the singular point is 3 (if £ 7^ 0) as can 
be checked using the invariants and tables given in 
I Artes et al, 2008| . 

In the normal form above, we consider the coef- 
ficient of the terms xy in both equations multiplied 
by 2 in order to make easier the calculations of the 
algebraic invariants we shall compute later. 

After an affine transformation and a time ho- 
mothety, a quadratic differential system with a 
semi-elemental triple node can be put into the fol- 
lowing normal form: 



X = 2xy + ky'^, 

y = y — x"^ + 2mxy + ny"^ 



where m, n and k are real parameters. We will 
study the bifurcations of system ^ in the 3- 
dimensional space with coordinates {m,n,k) G M^. 

Remark 2.2. After rescaling the time, we note that 
system ([7]) is symmetric in relation to the real pa- 
rameter k. So that, we will only consider k > 0. 



3. The Poincare compactification and the 
complex (real) foliation with singulari- 
ties on (MF2) 

A real planar polynomial vector field can be com- 
pactified on the sphere as follows. Consider the 
X, y plane as being the plane Z = 1 in the space 
with coordinates X, Y, Z. The central pro- 
jection of the vector field ^ on the sphere of ra- 
dius one yields a diffeomorphic vector field on the 
upper hemisphere and also another vector field on 
the lower hemisphere. There exists (for a proof see 
[Gonzales, 1969| ) an analytic vector field cp(^) on 
the whole sphere such that its restriction on the 
upper hemisphere has the same phase curves as the 
one constructed above from the polynomial vector 
field. The projection of the closed northern hemi- 
sphere H+ of §2 on Z = under {X,Y,Z) -)■ 
(X, Y) is called the Poincare disc. A singular point 
q of cp{^) is called an infinite (respectively, finite) 
singular point if g € S^, the equator (respectively, 
q S S^\S^). By the Poincare compactification of 
a polynomial vector field we mean the vector field 
cp{^) restricted to the upper hemisphere completed 
with the equator. 



Ideas in the remaining part of this section 
go back to Darboux's work [Darboux, 1878] . Let 
p{x, y) and q{x, y) be polynomials with real coeffi- 
cients. For the vector field 



d d 
P— + 



(8) 



dx ' dy' 

or equivalently for the differential system 

x=p{x,y), y = q{x,y), (9) 

we consider the associated differential 1-form 
01 = q{x, y)dx—p{x, y)dy, and the differential equa- 
tion 

01 = . (10) 
(7) Clearly, equation (jlOp defines a foliation with sin- 



gularities on C . The affine plane C is com- 
pactified on the complex projective space CP^ = 
(C^ \ {0}) / where {X, Y, Z) ~ {X' , Y\ Z') if and 
only if [X,Y,Z) = \(X' ,Y' , Z') for some complex 
A / 0. The equivalence class of {X, Y, Z) will be 
denoted hy[X:Y:Z]. 

The foliation with singularities defined by equa- 
tion (jlOp on can be extended to a foliation with 
singularities on CP^ and the 1-form 0i can be ex- 
tended to a meromorphic 1-form on CP^ which 
yields an equation = 0, i.e. 

A{X,Y, Z)dX + B{X,Y, Z)dY + C{X,Y, Z)dZ = 0, 

(11) 

whose coefficients A, B, C are homogeneous poly- 
nomials of the same degree and satisfy the relation: 

A{X, y, Z)X + B{X, Y, Z)Y + C(X, y, Z)Z = 0, 

(12) 

Indeed, consider the map i : \ {Z = 0} — )• C^, 
given by i{X, Y, Z) = (X/Z, Y/Z) = (x, y) and sup- 
pose that max{deg(p), deg(g)} = m > 0. Since 
x = X/Z and y = Y/Z we have: 

dx = {ZdX - XdZ) /Z^, dy = {ZdY - YdZ) /Z^, 

the pull-back form i*(0i) has poles at Z = and 
yields the equation 

r(0i) =q{X/Z,Y/Z){ZdX - XdZ)/Z'^ 

-p{X/Z,Y/Z){ZdY - YdZ)/Z'^ = 0. 

Then, the 1-form = Z™+2i*(0i) in \ {Z / 0} 
has homogeneous polynomial coefficients of degree 
m + 1, and for Z = the equations = and 
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i*(0i) = have the same solutions. Therefore, the 
differential equation = can be written as (jlip . 
where 

A{X, Y, Z) =ZQ{X, Y, Z) = Z'^+^iX/Z, Y/Z), 
B{X, Y,Z) = - ZP{X, Y, Z) = -Z'^+^p{X/Z, Y/Z), 
C{X, Y, Z) =YP{X, Y, Z) - XQ{X, Y, Z). 

(13) 

Clearly A, B and C are homogeneous polyno- 
mials of degree m + 1 satisfying (jl2p . 

In particular, for our quadratic systems ([7]), A, 
B and C take the following forms 

A{X, y, Z) =ZQ{X, Y, Z) = -Z{X'^ - 2mXY 

- nY"^ - YZ), 

B{X, Y,Z)=- ZP{X, Y, Z) = -YZ{2X + kY), 

CiX, Y, Z) =YP{X, Y, Z) - XQ{X, Y, Z) 

=X^ - 2mX^Y + 2XY^ (14) 

- nXY'^ + kY^ - XYZ. 

We note that the straight line Z = is always 
an algebraic invariant curve of this foliation and 
that its singular points are the solutions of the sys- 
tem: A{X,Y,Z) = B{X,Y,Z) = C{X,Y,Z) = 0. 
We note also that C{X,Y,Z) does not depend on 
b. 

To study the foliation with singularities defined 
by the differential equation (fTTj) subject to 
with A, B, C satisfying the above conditions in the 
neighborhood of the line Z = 0, we consider the 
two charts of CP^: {u,z) = {Y/X,Z/X), X / 0, 
and {v,w) = {X/Y, Z/Y), Y 0, covering this 
line. We note that in the intersection of the charts 
{x,y) = {X/Z,Y/Z) and {u,z) (respectively, {v,w)) 
we have the change of coordinates x = 1/z, y = u/z 
(respectively, x = v/w, y = 1/w). Except for the 
point [0 : 1 : 0] or the point [1:0:0], the foliation 
defined by equations (fTT]l . (fT2]) with A, B, C as in 
(jl3|) yields in the neighborhood of the line Z = 
the foliations associated with the systems 

tt =uP{l, u, z) — Q{l,u, z) = C(l, n, z), 
i =zP{l, n, z), 



(15) 



or 



i) =vQ{v, l,w) — P{v, l,w) = —C{v,l,w), 
w =wP{v, 1, w). 



(16) 



4. A few basic properties of quadratic sys- 
tems relevant for this study 

We list below results which play a role in the 
study of the global phase portraits of the real pla- 
nar quadratic systems ([T|) having a semi-elemental 
triple node. 

The following results hold for any quadratic 
system: 

(i) A straight line either has at most two (finite) 
contact points with a quadratic system (which 
include the singular points) , or it is formed by 
trajectories of the system; see Lemma 11.1 of 
|Ye et al, 1986| . We recall that by definition 
a contact point of a straight line L is a point of 
L where the vector field has the same direction 
as L, or it is zero. 

(ii) If a straight line passing through two real fi- 
nite singular points ri and r2 of a quadratic 
system is not formed by trajectories, then it is 
divided by these two singular points in three 
segments oorT, fir2 and r200 such that the 
trajectories cross ooFf and ¥200 in one direc- 
tion, and they cross fir2 in the opposite di- 



rection; see Lemma 11.4 of Ye et al., 1986 



In a similar way we can associate a real foliation 
with singularities on MP^ to a real planar polyno- 
mial vector field. 



(iii) If a quadratic system has a limit cycle, then 
it surrounds a unique singular point, and this 
point is a focus; see [Coppel, 1966| . 

Theorem 4.1. Any graphic or degenerated graphic 
in a real planar polynomial differential system must 
either 

1 ) surround a singular point of index greater than 
or equal to +1, or 

2) contain a singular point having an elliptic sector 
situated in the region delimited by the graphic, or 

3) contain an infinite number of singular points. 

Proof. Let 5 be a simply connected invariant 
closed set under the fiow of a vector field. In 
[Artes et al., 1998] the index of dS is given by: 
Yll^=\{Ei - Hi + l)/2, where Ei (respectively. Hi) 
is the number of elliptic (respectively, hyperbolic) 
sectors which are inside the region delimited by S 
of the singular points forming the border. Also the 
index of S is given by the index of dS plus the sum 
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of the indices of the singular points in the interior 
of S. 

Prom the same paper, Proposition 4.8 claims 
that given a vector field X or p{X) and S an in- 
variant region topologically equivalent to (the 
closed disk) containing a finite number of singular 
points (both in dS or its interior), then the index 
of S is always +1. 

Now, assume that we have a graphic of a poly- 
nomial system. If it contains an infinite number 
of singular points (either finite or infinite) we are 
done. Otherwise, such a graphic altogether with 
its interior is an invariant region as defined in 
[Artes et al, 19981 must have index Since 
the index is positive, we must have some element, 
either in the interior or on the border which makes 
the index positive, and this implies the existence of 
either a point of index greater than or equal to 
in its interior or at least one elliptic sector coming 
from a singular point on the border and situated in 
the region delimited by the graphic. ■ 

5. Some algebraic and geometric concepts 

In this article we use the concept of intersection 
number for curves (see [Pulton, 1969| ). Por a quick 
summary see Sec. 5 of [Artes et al., 2006] . 

We shall also use the concepts of zero- 
cycle and divisor (see [Hartshorne, 1977] ) 
as specified for quadratic vector fields in 
[Schlomiuk et al., 2001 [ . Por a quick summary see 
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We shall also use the concepts of algebraic in- 
variant and T-comitant as used by the Sibirskii 
school for differential equations. Por a quick sum- 
mary see Sec. 7 of [Artes et al., 2006] . 

The invariants which are relevant in this study 
are among the ones which were relevant in the study 
of QW2 and will be described in the next section. 

6. The bifurcation diagram of the systems 
with a semi elemental triple node 

6.1. Bifurcation surfaces due to the 
changes in the nature of singularities 

Por systems d?]) we will always have (0, 0) as a finite 
singular point, a semi-elemental triple node. 

Prom Sec. 7 of [Artes et al, 2008[ we get the 
formulas which give the bifurcation surfaces of sin- 



gularities in M^^, produced by changes that may oc- 
cur in the local nature of finite singularities. Prom 
IS chlomiuk et al., 2005[ we get equivalent formulas 
for the infinite singular points. These bifurcation 
surfaces are all algebraic and they are the follow- 
ing: 

Bifurcation surfaces in M'^ due to multiplici- 
ties of singularities 

This is the bifurcation surface due to mul- 
tiplicity of infinite singularities as detected by 
the coefficients of the divisor D^{P,Q;Z) = 
Y.w&{z=o}nC¥^^w{P,QW , (here Iw{P,Q) de- 
notes the intersection multiplicity of P = with 
Q = at the point W situated on the line at infin- 
ity, i.e. Z = 0) whenever deg{{D^{P,Q; Z))) > 0. 
This occurs when at least one finite singular point 
collides with at least one infinite point. More 
precisely this happens whenever the homogenous 
polynomials of degree two, p2 and (72 in p and q 
have a common root. In other words whenever 
12 = ReSxip2, Q2)/y^ = 0. The equation of this sur- 
face is 

pL = k"^ + 4km — 4n = 0. 

(55)0 This is the bifurcation surface due to 
multiplicity of infinite singularities as de- 
tected by the coefficients of Dc{C,Z) = 
Ylw&{z=o}nC¥2lw{C,Z)W, i.e. this bifurca- 
tion occurs whenever at a point W of intersection 
of C = with Z = we have Iw{C, Z) > 2, i.e. 
when at least two infinite singular points collide 
at W . This occurs whenever the discriminant of 
C2 = C7(X,y,0) = Yp2{X,Y) - Xq2{X,Y) is zero 
where by p2 , q2 we denoted the second degree terms 
in p, q. We denote by r/ this discriminant. The 
equation of this surface is 

7? = - 32 - 27A;^ - 72fcm + IQrr? + 32km^ + A8n+ 
36kmn — 16m n — 24n -|- 4m n + 4n =0. 

C°° bifurcation surface in due to a strong 
saddle or a strong focus changing the sign of 
their traces (weak saddle or weak focus) 



^The numbers attached to these bifurcations surfaces do 
not appear here in increasing order. We just kept the same 
enumeration used in [Artes et al, 2006] to maintain coher- 
ence even though some of the numbers in that enumeration 
do not occur here. 
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{S3) This is the bifurcation surface due to the weak- 
ness of finite singularities, which occurs when the 
trace of a finite singular point is zero. The equation 
of this surface is given by 

% = 8 + k^ + 4n = 0, 

where 71 is defined in [Vulpe, 2011 j . This Ti is an 
invariant. 

C°° bifurcation surface in due to a node 
becoming a focus 

(56 ) This surface will contain the points of the pa- 
rameter space where a finite node of the system 
turns into a focus. This surface is a C°° but not 
a topological bifurcation surface. In fact, when 
we only cross the surface (Sq) in the bifurcation 
diagram, the topological phase portraits do not 
change. However, this surface is relevant for iso- 
lating the regions where a limit cycle surrounding 
an antisaddle (different from the triple node) can- 
not exist. Using the results of [Artes et ai, 2008] , 
the equation of this surface is given by W4 = 0, 
where 

1^4 = 64 48A;^ + k^ + 128km - 64n + 8k^n + 16n^. 

These are all the bifurcation surfaces of singu- 
larities of the systems ([7]) in the parameter space 
and they are all algebraic. We shall discover an- 
other bifurcation surface not necessarily algebraic 
and on which the systems have global connection 
of separatrices. The equation of this bifurcation 
surface can only be determined approximately by 
means of numerical tools. Using arguments of con- 
tinuity in the phase portraits we can prove the ex- 
istence of this not necessarily algebraic component 
in the region where it appears, and we can check it 
numerically. We will name it the surface (^7). 

Remark 6.1. Even though we can draw a 3- 
dimensional picture of the algebraic bifurcation sur- 
faces of singularities in (see Fig. [2] for an ex- 
ample), it is pointless to try to see a single 3~ 
dimensional image of all these four bifurcation sur- 
faces together in the space M^. As we shall see later, 
the full partition of the parameter space obtained 
from all these bifurcation surfaces has 63 parts. 

Due to the above remark we shall foliate the 
3-dimensional bifurcation diagram in by planes 
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Fig. 2. The 3-dimensional picture of the surface [Sq) 
(when a finite node becomes a focus). 

k = ko, kf) constant. We shall give pictures of the 
resulting bifurcation diagram on these planar sec- 
tions on an affine chart on M^. In order to detect 
the key values for this foliation, we must find the 
values of parameters where the surfaces intersect to 
each other. As we mentioned before, we will be only 
interested in non-negative values of k to construct 
the bifurcation diagram. 

The following set of seven results study the sin- 
gularities of each surface and the simultaneous in- 
tersection points of the bifurcation surfaces, or the 
points or curves where two bifurcation surfaces are 
tangent. 

As the final bifurcation diagram is quite com- 
plex, it is useful to introduce colors which will be 
used to talk about the bifurcation points: 

(a) the curve obtained from the surface (Si) is 
drawn in blue (a finite singular point collides 
with an infinite one); 

(b) the curve obtained from the surface (^3) is 
drawn in green (when the trace of a singular 
point becomes zero); 

(c) the curve obtained from the surface (^5) is 
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drawn in red (two infinite singular points col- 
lide); 

(d) the curve obtained from the surface (Sq) is 
drawn in black (an antisaddle different from the 
origin is on the verge of turning from a node to 
a focus or vice versa); and 

(e) the curve obtained from the surface (Sj) is 
drawn in purple (the connection of separatri- 
ces). 

Lemma 6.2. Concerning the singularities of the 
surfaces, it follows that: 

(i) (Si) and (S3) have no singularities; 

(a) (S^) has a curve of singularities given by 
4m? + 3n - 6 = 0; 

(Hi) (Sq) has a singularity on the straight line 
(m, 2, 0) on slice k = 0. Besides, this sur- 
face restricted to k = is part of the surface 

(s,). 

Proof. It is easy to see that the gradient of (»Si) 
and (iSs) is never null for all {m,n,k) G M^; so (i) 
is proved. In order to prove (ii) we compute the 
gradient of rj and we verify that it is null whenever 
m = -3\^/2 and n = 2 - 3v^, for all A; > 0. It 
is easy to see that these values of m and n for all 
A; > lie on the curve 4m^ + 3n — 6 = 0. Finally, 
considering the gradient of the surface (Sq), it is 
identically zero at the point (0, 2, 0) which lies on 
the straight line (m, 2, 0) whenever k = 0. More- 
over, if /c = 0, we see that the equation of (Sq) is 
(—2 + n)^, which is part of (^5), proving (Hi). ■ 

Lemma 6.3. Ifk = 0, the surfaces (Si) and (S3) 
have no intersection. For all k > 0, they intersect 
in the point {-{A + k^)/2k, -2 - k'^/A, k). 

Proof. By solving simultaneously both equations of 
the surfaces [Si) and (^3) for all A; > 0, we obtain 
the point {-{A+k'^)/2k,-2-k'^ /A^k). We also note 
that if A = there is no intersection point. ■ 

Lemma 6.4. If k = 0, the surfaces (Si) and (S^) 
intersect at the points (— \/2, 0, 0) and (\/2, 0, 0). 
For all k > 0, they intersect along the curve 
71 (m,n) = —64 + 32m^ + 16n — = and they 



have a 2-order contact along the curve 72(m,n) = 
1 + 2m2 + 2n + rn^n + = 0. 

Proof. By solving simultaneously both equations 
of the surfaces {Si) and (^5) for fc = 0, we ob- 
tain the two solutions mi = — \/2, ni = and 
m,2 = \f2, 71,2 = 0, proving the first part of the 
lemma. For all fc > 0, the simultaneous solu- 
tions of the equations are the three points: ri = 
(-\/2-A/4,-v^A:, A;), = (\/2 - A/4, \/2A, A) and 
r3 = (-(4 + A^)/2A, -2 - A^/4, A). By computing 
the resultant with respect to A of {Si) and (5.5), we 
see that Resfc[(5i), (55)] = — 167i(m, n)(72(m, /i))^, 
where 71 (m,n) and 72 (m,n) are as stated in the 
statement of the lemma. It is easy to see that 
71 (m,n) has two simple roots which are ri and r2, 
and ra is a double root of (72(777., n))^. So that, 
the surfaces intersect transversally along the curve 
71 (m, n) and they have a 2-order contact along the 
curve 72 (m, ra). ■ 

Lemma 6.5. If k = 0, the surfaces (Si) and (Sq) 
have no intersection. For all A > 0, they have a 2- 
order contact along the curve l + 2m^ + 2n + m^n-|- 
n2 = 0. 

Proof. By solving the system formed by the equa- 
tions of the surfaces {Si) and {Sq), we find the point 
r = (-(4+A;2)/2A, -2-A:V4, A), for all A > 0, which 
lies on the curve 1 + 2?7i^ + 2n + m?n + = 0. Wc 
claim that the surfaces {Si) and {Sq) have a 2-ordcr 
contact point at r. Indeed, we have just shown that 
the point r is a common point of both surfaces. Ap- 
plying the change of coordinates given by n = (v -|- 
Am+A^)/4, G R, we see that the gradient vector of 
{Si) is V/i(r) = (0,0,0) while the gradient vector 
of (56) is VW4(r) = (0,0,8(-4 + 16/P + SA;^)), 
whose last coordinate is always positive for all 
A > 0. As it does not change its sign, the vec- 
tor VW4^{r) will always point upwards in relation 
to {Si) restricted to the previous change of coor- 
dinates. Then, the surface {Sq) remains only on 
one of the half-spaces delimited by the surface {Si), 
proving our claim. ■ 

Lemma 6.6. If k = 0, the surfaces (S3) and (S^) 
intersect at the points (—2,-2,0) and (2,-2,0). 
For all k > 0, they intersect at the points ri = 
{{32k -k^- 7(64- A2)3) /256, -2 - k'^/A, k),r2 = 
(-(4 + k'^)/2k, -2 - k) and = ((32A - k^ + 
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v/(64 - A;2)3)/256, -2 - k'^/A, k). 

Proof. The result follows easily by solving the sys- 
tem formed by the equations of the surfaces. ■ 

Corollary 6.7. If k = 2^/2, the points ri and r2 
of Lemma \6.6\ are equal and they correspond to the 
singularity of the surface (S^). 

Proof. Replacing k = 2\f2 at the expressions of the 
points ri, r2 and r3 described in Lemma [6 .Gt we see 
that r\ = r2 and they are equal to the singularity 
(-3\/2/2, -4, 2^/2) of the surface {S^). ■ 

Remark 6.8. We observe that the values k = Q and 
k = 2^/2 win be very important to describe the 
bifurcation diagram due to the "rich" change on 
the behavior of the surfaces. 

Lemma 6.9. If k = 0, the surfaces (S^) and (Sq) 
intersect along the straight line (m, 2,0), for all 
m G M. For all k > 0, they have a 2-order con- 
tact point at (-(4 + k'^)/2k, -2 - P/4, k). 

Proof. Replacing = in the equations of the 
surfaces and solving them in the variables m and 
n, we find that m G M and n = 2, implying 
the existence of intersection along the straight line 
(m, 2,0), m E M. For all k > 0, the solution 
of the equations of the surfaces is the point r = 
(-(4 + k'^)/2k,-2 - k'^/A,k). We claim that the 
surfaces (^5) and (Sq) have a 2-order contact point 
at r. We shall prove this claim by showing that 
each one of the surfaces (^5) and (Sq) remains on 
only one of the half-spaces delimited by the plane 
(5i) and their unique common point is r. Indeed, it 
is easy to see that the point r is a common point of 
the three surfaces. By applying the change of coor- 
dinates given hy n = (v + km + A;^)/4, t; G R, as in 
the proof of Lemma [63t we see that the surface (Sq) 
remains on only one of the half-spaces delimited by 
the plane (Si). On the other hand, numerical calcu- 
lations show us that the surface (^5) is zero valued 
around the point r and it assumes negative values 
otherwise, showing that {S5) remains on the other 
half-space delimited by the plane (5i). ■ 

Lemma 6.10. The curve r{k) = {—3\/k/2, 
2 — 3-Vk^,k) of (S5) (i.e., its set of singulari- 
ties) cannot belong to the region where W4 > 



and fj, < 0. 

Proof. We consider the real continuous function 
g = (^H^4)L,(,) = (^-2)' (^ + 6)^, 

whose zeroes are and 2V2. It is easy to see that 
g is always positive in {0,2\/2) U (2\/2, 00), imply- 
ing that the functions fj, and W4 calculated at r{k) 
cannot have different signs, proving the lemma. ■ 

Now we shall study the bifurcation diagram 
having as reference the values of k where significant 
phenomena occur in the behavior of the bifurcation 
surfaces. 

According to the Remark 16. 8| these values are 
k = and k = 2\f2. So, we only need to add two 
more slices with some intermediate values. 

We take then the values: 

fco = 0, A;i = 1, 
k2 = 2V2, k-i = 3. 

The values indexed by positive even indices cor- 
respond to explicit values of k for which there is a 
bifurcation in the behavior of the systems on the 
slices. Those indexed by odd ones are just interme- 
diate points (see Figs. [3] to ED . 

We now describe the labels used for each part. 
The subsets of dimensions 3,2,1 and 0, of the parti- 
tion of the parameter space will be denoted respec- 
tively by V , S, L and P for Volume, Surface, Line 
and Point, respectively. The surfaces are named 
using a number which corresponds to each bifurca- 
tion surface which is placed on the left side of S*. To 
describe the portion of the surface we place an in- 
dex. The curves that are intersection of surfaces are 
named by using their corresponding numbers on the 
left side of L, separated by a point. To describe the 
segment of the curve we place an index. Volumes 
and Points are simply indexed (since three or more 
surfaces may be involved in such an intersection). 

We consider an example: the surface {Si) splits 
into 5 different two-dimensional parts labeled from 
l^i to 155, plus some one-dimensional arcs labeled 
as l.iLj (where i denotes the other surface inter- 
sected by {Si) and j is a number), and some zero- 
dimensional parts. In order to simplify the labels 
in Figs. [To] to [13] we see VI which stands for the 
notation Vi. Analogously, ISl (respectively, 
1.2L1) stands for ISi (respectively, 1.2Li). And 
the same happens with many other pictures. 



Global phase portraits of quadratic polynomial differential systems with a semi-elemental triple node 11 



Some bifurcation surfaces intersect on = 
or have singularities there. The restrictions of the 
surfaces on /c = are: the surface (S^) has a sin- 
gularity at the point (0, 2, 0) and it is the union of 
a parabola and a straight line of multiplicity two, 
which in turn coincides with the bifurcation surface 
(Sq); the surface coincides with the horizon- 
tal axis and the bifurcation surface (S^) becomes a 
straight line parallel to the horizontal line having 
intersection points only with the surface (S^). 

As an exact drawing of the curves pro- 
duced by intersecting the surfaces with slices 
gives us very small regions which are difficult 
to distinguish, and points of tangency are al- 
most impossible to recognize, we have produced 
topologically equivalent figures where regions are 
enlarged and tangencies are easy to observe. The 
reader may find the exact pictures in the web page 
http: / / mat.uab.es/~artes/articles/qvftn/qvftn.html 

As we increase the value of k, other changes in 
the bifurcation diagram happen. When k = 1, the 
surface {S^) has two connected components and a 
cusp point as a singularity which remains on the 
left side of the surface (Sq) until k = 2\/2 (see Fig. 
H]). At this value, the cusp point is the point of 
contact among all the surfaces, as we can see in 
Fig. El and when A; = 3, the cusp point of (^5) 
lies on the right side of surfaces (Sq) and (Si) (see 
Fig. O and Lemma l6.1Up . In order to comprehend 
that the "movement" of the cusp point of the sur- 
face (55) implies changes that occur in the bifurca- 
tion diagram, we see that when k = 1 we have a 
"curved triangular" region formed by the surfaces 
(53) and (55), the cusp point of (^5) and the points 
of intersection between both surfaces. The "trian- 
gle" bounded by these elements yields 15 subsets: 
three 3-dimensional subsets, seven 2-dimensional 
ones and five l~dimensional ones. In Fig. O we 
see that the "triangle" has disappeared and it has 
become a unique point which corresponds to the 
point of contact of all the surfaces and the cusp 
point of surface {S5). Finally, when = 3, the "tri- 
angle" reappears and yields also 15 subsets of same 
dimensions, but different from the previous ones. 

All other regions of the parameter space related 
to singular points remain topologically the same 
with respect to the algebraic bifurcations of singu- 
larities when moving from Figs. [UtoEl 

We recall that the black curve (Sq) (or W4) 
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Fig. 3. Slice of parameter space when k = 0. 




Fig. 4. Slice of paranieter space when fc = 1. 
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Fig. 5. Slice of parameter space when k = 2^/2. 
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Fig. 6. Slice of parameter space when A; = 3. 



means the turning of a finite antisaddle different 
from the origin from a node to a focus. Then, ac- 
cording to the general results about quadratic sys- 
tems, we could have limit cycles around such focus 
for any set of parameters having VF4 < 0. 

Remark 6.11. Wherever two parts of equal dimen- 
sion d are separated only by a part of dimension 
d — 1 of the black bifurcation surface (Sq), their 
respective phase portraits are topologically equiva- 
lent since the only difference between them is that 
a finite antisaddle has turned into a focus without 
change of stability and without appearance of limit 
cycles. We denote such parts with different labels, 
but we do not give specific phase portraits in pic- 
tures attached to Theorem 11.11 for the parts with 
the focus. We only give portraits for the parts with 
nodes, except in the case of existence of a limit cycle 
or a graphic where the singular point inside them is 
portrayed as a focus. Neither do we give specific in- 
variant description in Sec. [7] distinguishing between 
these nodes and foci. 

6.2. Bifurcation surfaces due to connec- 
tions 

We now place for each set of the partition on = 3 
the local behavior of the flow around all the sin- 
gular points. For a specific value of parameters of 
each one of the sets in this partition we compute the 
global phase portrait with the numerical program 
P4 [Dumortier et al., 2006 j . In fact, many (but not 
all) of the phase portraits in this work can be ob- 
tained not only numerically but also by means of 
perturbations of the systems of codimension one. 

In this slice we have a partition in 2- 
dimensional regions bordered by curved polygons, 
some of them bounded, others bordered by infinity. 
From now on, we use lower letters provisionally to 
describe the sets found algebraically so not to inter- 
fere with the final partition described with capital 
letters. For each 2~dimensional region we obtain 
a phase portrait which is coherent with those of 
all their borders. Except one region. Consider the 
segment 3si in Fig. [6l On it we have a weak fo- 
cus and a Hopf bifurcation. This means that either 
in ^5 or vq we must have a limit cycle. In fact it 
is in vq. The same happens in 3s2, so a limit cy- 
cle must exist either in V14 or vj. However, when 
approaching 6si or 6s2, this limit cycle must have 
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disappeared. So, either f 7 or f 14 must be splitted in 
two regions separated by a new surface (S^) having 
at least one element 7Si such that one region has 
limit cycle and the other does not, and the border 
7Si must correspond to a connection between sepa- 
ratrices. Numerically it can be checked that it is the 
region vj the one which splits in V7 without limit 
cycles and V15 with one limit cycle. It can also be 
analytically proved (see Proposition 16.12]) that the 
segment 5s4 must be splitted in two segments 5S4 
and by the 1-dimensional subset 5.7 Li. The 
other border of 75i must be 1.3Li for coherence. 
We plot the complete bifurcation diagram in Fig. 
[T3I We also show the sequence of phase portraits 
along these subsets in Fig. [71 

Notice that the limit cycle which is "born" by 
Hopf on 3Si either "dies" on 5S'4 or "survives" when 
crossing 7/ = 0, if we do it through bS^, and then it 
"dies" either on 7 Si or again by Hopf in 3^2. 

Surface {S7), for a concrete k > 2\/2, is a curve 
which starts on l.SLi and may either cut 5s4, or 
not. We are going to prove that, at least for a 
concrete fco, (Si) must cut it, and consequently it 
must do the same for an open interval around ko, 
thus proving the existence of subsets 55*4, 5^5 and 
5.7Li which have different phase portraits. 

Proposition 6.12. The following statements hold: 

(a) System 1^ with {m,n,k) = (-29/2,-105/4,7) 
has an even number of limit cycles (counting 
their multiplicities), and possibly this number 
is zero; 

(h) System ([7]) with {m,n,k) = 
(-49/2,-185/4,12) has an odd number of 
limit cycles (counting their multiplicities), and 
possibly this number is one; 

Proof, (a) We see that the system with rational 
coefficients 

X = 2xy + 7y'^, , . 

y = y _ _ 29xy - i25y2 U*) 

is a representative of the red surface (^5) which 
belongs to the subset 5s4. 

We have to show that there exists a hyperbola 

7i = ax^ + bxy + y'^ + dx + ey + f = 

which isolates the focus of (jlSp on the region where 
H < and x > 0, and with the property that at 




Fig. 7. Sequence of phase portraits in slice fc = 3. We 
start from V5. We recall that the phase portrait 35*1 
is equivalent to the phase portrait V5 up to a weak fo- 
cus (represented by a little black square) in place of the 
focus. When crossing 3si, we shall obtain the phase 
portrait Vg in subset vq . From this point we may choose 
three different ways to reach the subset by crossing 
5S4: (1) from the phase portrait 55*4 to the P7; (2) from 
the phase portrait 5.7Li to the 751; and (3) from the 
phase portrait 53^ to the V15, from where we can move 

to V14. 
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Fig. 8. The hyperbola m the phase portrait ol 5^4 



Fig. 9. The hyperbola in the phase portrait of bS^. 



each of its points the flow crosses the hyperbola in 
only one direction, as we can see in Fig. [8j By 
proving the existence of this hyperbola, we shall 
prove that (jlSp has an even number of limit cycles. 

For convenience and making easier the calcula- 
tions, we impose that the hyperbola passes through 
two infinite singular points of (jlSp with the same 
tangencies of the affine separatrices. With all these 
features we have just one free parameter which is 
used to force the hyperbola to pass through a con- 
crete finite point. In resume, we take a = 1/14, 
b = 57/28, 



d 



1 

392 

/ 



148225e2 + 1 + 399e - 1 



e2 - 28cle + {2M - ef 



1516A/(28d - e)2 - 110 

and e = -324/10000. 

This hyperbola has a component fully included 
in the fourth quadrant and it is easy to check that 
the scalar product of its tangent vector with the 
flow of the vector field does not change its sign and 
the flow moves outwards the region Ti < 0. Since 
the focus is repellor, this is consistent with the ab- 
sence of limit cycles (or with an even number of 
them, counting their multiplicities). 

(b) We see that the system with rational coef- 
ficients 



2xy + I2y\ 
y — — 49xy 



185„,2 

4 y 



(19) 



is a representative of the red surface {S^) which 
belongs to the subset 5s4. 

Analogously, we have to show that there exists 
a hyperbola 

% = ax^ + bxy + y"^ + dx + ey + f = 



which isolates the focus of ()19p on the region where 
H < and x > 0, and with the property that at 
each of its points the flow crosses the hyperbola in 
only one direction, as we can see in Fig. [9l By 
proving the existence of this hyperbola, we shall 
prove that (fT9]l has an odd number of limit cycles. 

By using the same technique as before, we take 
a = 1/24, 6 = 97/48, 

1 



1152 



\/l299600e2 + 1 + 1164e - 1 



/ 



- 48de + ev^(48d - e)2 



4516^(48d - e)2 - 190 

and e = -18663/100000. 

This hyperbola has a component fully included 
in the fourth quadrant and it is easy to check that 
the scalar product of its tangent vector with the 
flow of the vector field does not change its sign and 
the flow moves inwards the region T-L < 0. Since the 
focus is repellor, this is consistent with the presence 
of one limit cycle (or with an odd number of them, 
counting their multiplicities). ■ 

We cannot be sure that this is all the addi- 
tional bifurcation curves in this slice. There could 
exist others which are closed curves which are small 
enough to escape our numerical research. For all 
other two-dimensional parts of the partition of this 
slice whenever we join two points which are close 
to two different borders of the part, the two phase 
portraits are topologically equivalent. So we do not 
encounter more situations than the one mentioned 
above. 

As we vary k in (2\/2, oo), the numerical re- 
search shows us the existence of the phenomenon 
just described, but for the values of k in [0,2\/2), 
we have not found the same behavior. 
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Fig. 10. Complete bifurcation diagram for slice k ~ 0. 

In Figs. [To] to [13] we show the complete bifurca- 
tion diagrams. In Sec. [7] the reader can look for the 
topological equivalences among the phase portraits 
appearing in the various parts and the selected no- 
tation for their representatives in Fig. [TJ 

6.3. Other relevant facts about the bifurca- 
tion diagram 

The bifurcation diagram we have obtained for QTN 
is completely coherent. By this, we mean that if we 
take any two points in the parameter space and join 
them by a continuous curve, along this curve the 
changes in phase portraits that occur when crossing 
the different bifurcation surfaces we mention can be 
completely explained. 

However, we cannot be sure that this bifurca- 
tion diagram is the complete bifurcation diagram 
for QTN due to the possibility of "islands" inside 
the parts bordered by unmentioned bifurcation sur- 
faces. In case they exist, these "islands" would not 
mean any modification of the nature of the singu- 
lar points. So, on the border of these "islands" we 
could only have bifurcations due to saddle connec- 
tions or multiple limit cycles. 
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Fig. 11. Complete bifurcation diagram for slice k = 1. 




Fig. 12. Complete bifurcation diagram for slice k — 
2v/2. 
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Fig. 13. Complete bifurcation diagram for slice k — 3. 



In case there were more bifurcation surfaces, we 
should still be able to join two representatives of 
any two parts of the 63 parts found until now with 
a continuous curve either without crossing such bi- 
furcation surface or, in case the curve crosses it, it 
must do it an even number of times without tan- 
gencies, otherwise one must take into account the 
multiplicity of the tangency, so the total number 
must be even. This is why we call these potential 
bifurcation surfaces ^Hslands" . 

To give an example of such a potential "island" , 
we consider region Vi where we have a phase por- 
trait having a finite antisaddle, a saddle and two 
pairs of infinite antisaddles and one pair of infinite 
saddles. This phase portrait is topologically equiva- 
lent (modulo limit cycles and taking the triple node 
as a simple antisaddle) with the phase portrait 9.1 
from [Artes et al., 1998] where all structurally sta- 
ble quadratic vector fields were studied, (see the 
first phase portrait of Fig. [T^ . 

We note that in [Artes et al., 1998] it is proved 
that structurally stable (modulo limit cycles) 
quadratic vector fields can have exactly 44 differ- 
ent phase portraits. In the case of system ([7]), we 
have a semi-elemental triple node which topologi- 
cally behaves like an elemental node, and the phase 



Fig. 14. Example of a potential "island" . 

portraits in generic regions on bifurcation diagram 
will look like structurally stable ones. From those 
44, two have no singular points, one has no finite 
antisaddles and 33 have four finite singular points, 
so obviously they cannot appear in QTN. From 
the remaining 8, only 7 appear in our description 
of QTN. There are two potential reasons for the 
absence of the remaining case: (1) it cannot be re- 
alized within QTN, or (2) it may live in such "is- 
lands" where the conditions for the singular points 
are met, but the separatrix configuration is not the 
one that we have detected as needed for the coher- 
ence. 

For example, the structurally stable phase por- 
trait 9.2 has so far not appeared anywhere, but it 
could perfectly fit in an "island" inside Vi (or Vu) 
where we have phase portrait 9.1. The transition 
from 9.1 to 9.2 consists on the existence of a het- 
eroclinic connection between the finite saddle and 
one of the infinite saddles as it can be seen in Fig. 
[Hi We also show (in the middle of this figure) the 
unstable phase portrait from which could bifurcate 
and also has the potential to be on the bifurcation 
surface delimiting the "island" . 

7. Completion of the proof of the main the- 
orem 

In the bifurcation diagram we may have topolog- 
ically equivalent phase portraits belonging to dis- 
tinct parts of the parameter space. As here we 
have 63 distinct parts of the parameter space, 
to help us identify or to distinguish phase por- 
traits, we need to introduce some invariants and 
we actually choose integer-valued invariants. All 
of them were already used in [Llibre et al., 2004 



Artes et al, 2006] . These integer-valued invariants 
yield a classification which is easier to grasp. 

Definition 7.1. We denote by Ii{S) the number 
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of the real finite singular points. This invariant is 
also denoted by Nirj(S') |Artes et al, 2006| . 

Definition 7.2. We denote by l2{S) the sum 
of the indices of the real finite singular points. 
This invariant is also denoted by deg{DIf{S)) 
Artes et al, 2006| . 

Definition 7.3. We denote by I^iS) the number 
of the real infinite singular points. This invariant is 
also denoted by Nk^oo(5') [Artes et al, 2006] . 

Definition 7.4. We denote by l4{S) the sequence 
of digits (each one ranging from to 4) such that 
each digit describes the total number of global or 
local separatrices (different from the line of infinity) 
ending (or starting) at an infinite singular point. 
The number of digits in the sequences is 2, 4 or 6 
according to the number of infinite singular points. 
We can start the sequence at anyone of the infinite 
singular points but all sequences must be listed in 
a same specific order either clockwise or counter- 
clockwise along the line of infinity. 

In our case we have used the clockwise sense 
beginning from the top-most infinite singular point 
in the pictures shown in Fig. [TJ 

Definition 7.5. We denote by h{S) a digit which 
gives the number of limit cycles. 

As we have noted previously in Remark 16. IH 
we do not distinguish between phase portraits 
whose only difference is that in one we have a fi- 
nite node and in the other a focus. Both phase 
portraits are topologically equivalent and they can 
only be distinguished within the class. In case 
we may want to distinguish between them, a new 
invariant may easily be introduced. 

Theorem 7.6. Consider the family QTN of all 
quadratic systems with a semi-elemental triple 
node. Consider now all the phase portraits that we 
have obtained for this family. The values of the 
affine invariant I = (Ii, I2, 13, 14, 15) given in the 
following diagram yield a partition of these phase 
portraits of the family QTN . 

Furthermore, for each value of I in this dia- 
gram there corresponds a single phase portrait; i.e. 
S and S' are such that I{S) = I{S'), if and only if 



S and S' are topologically equivalent. 

The bifurcation diagram for QTN has 63 parts 
which produce 28 topologically different phase por- 
traits as described in Table 1. The remaining 35 
parts do not produce any new phase portrait which 
was not included in the 28 previous. The difference 
is basically the presence of a strong focus instead of 
a node and vice versa. 

The phase portraits having neither limit cy- 
cle nor graphic have been denoted surrounded by 
parenthesis, for example (Vi); the phase portraits 
having one limit cycle have been denoted sur- 
rounded by brackets, for example [Vq]; the phase 
portraits having one graphic have been denoted sur- 
rounded by {}, for example {5 5*4}. 

Proof. The above result follows from the results 
in the previous sections and a careful analysis of 
the bifurcation diagrams given in Sec. [6] in Figs. 
[TOl to \T3\ the definition of the invariants Ij and 
their explicit values for the corresponding phase 
portraits. ■ 

We first make some observations regarding the 
equivalence relations used in this work: the affine 
and time rescaling, and topological equivalences. 

The coarsest one among these three is the topo- 
logical equivalence and the finest is the affine equiv- 
alence. We can have two systems which are topo- 
logically equivalent but not -equivalent. For ex- 
ample, we could have a system with a finite anti- 
saddle which is a structurally stable node and in 
another system with a focus, the two systems being 
topologically equivalent but belonging to distinct 
-equivalence classes, separated by a surface {Sq 
in this case) on which the node turns into a focus. 

In Table 2 we listed in the first column 28 parts 
with all the distinct phase portraits of Fig. [TJ Cor- 
responding to each part listed in column 1 we have 
in its horizontal block, all parts whose phase por- 
traits are topologically equivalent to the phase por- 
trait appearing in column 1 of the same horizontal 
block. 

In the second column we have put all the parts 
whose systems yield topologically equivalent phase 
portraits to those in the first column but which may 
have some algebro-geometric features related to the 
position of the orbits. 
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In the third (respectively, fourth, and fifth) col- 
umn we list all parts whose phase portraits have 
another antisaddle which is a focus (respectively, a 
node which is at a bifurcation point producing foci 
close to the node in perturbations, a node-focus to 
shorten, and a finite weak singular point). 

Whenever phase portraits appear on a horizon- 
tal block in a specific column, the listing is done 
according to the decreasing dimension of the parts 
where they appear, always placing the lower dimen- 
sions on lower lines. 
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Table 1: Geometric classification. 
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Table 2: Topological equivalences. 
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Fig. 15. The perturbations of phase portrait Vu yield- 
ing the structurally unstable phase portraits A22 and 

^23- 

7. 1 . Proof of the main theorem 

The bifurcation diagram described in Sec. [6l plus 
Table 1 of the geometrical invariants distinguish- 
ing the 28 phase portraits, plus Table 2 giving 
the equivalences with the remaining phase portraits 
lead to the proof of the main statement of Theorem 

o ^ 

In [Artes fc Llibre, 2013| the authors are study- 
ing all phase portraits of quadratic systems hav- 
ing exactly one saddle-node or one connection of 
separatrices. By using a similar technique as the 
one used in [Artes et al, 1998] for the structurally 
stable ones, they have produced a complete list of 
topologically possible structurally unstable systems 
of codimension one (modulo limit cycles), they have 
erased many of them proving their impossibility and 
they have proved the existence of many others (180 
just before this paper), and it remains 24 which 
escaped up to now both the proof of their impossi- 
bility and finding an example. 

Our system in QTN Vu yields an example of 
their "wanted" case A23 by perturbating the triple 
node while producing the desired phase portrait as 
may be seen in Fig. [T5l 
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